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OF DAMPING IN FLUTTER MODELS 
By Robert P.- Coleman 

SUMMARY 



The problem of determining values of structural damp- 
ing for use in flutter calculations -is discussed. The. 
concept of equivalent -viscous damping is reviewed and its 
relation to the structural damping coefficient g intro- 
duced in N.A.C.A- Technical Report No. 685 is shown. The 
theory of normal modes is reviewed and a number of methods 
are described for separating the motions associated with 1 
different modes. Equations are developed for use in eval- 
uating the damping parameters from experimental data. 

Experimental results of measurements of . damping in 
several flutter models are presented. 



INTRODUCTION 



One important step. in the study of the flutter prop- 
erties of an airplane structure is the determination of 
the' structural damping. ■ In an investigation of flutter 
carried out in the N.A.Q.A. 8-foot hi?h-speed tunnel, it ■ 
was- desirable to determine the damping in the models test- 
ed. The present report is a description of the methods 
used, together with a review and a critical discussion of 
the principles and::the derivations pertaining to the meas- 
urement of damping parameters.- 

The presence of damping in a structure can be in- 
ferred from and its amount can be measured by a number of 
different effects.- The principal effects that depend upon 
damping are : 

(1) The rate of decay of free vibrations'. 

(2) The amplitudes produced by given applied forces 

at a resonant frequency. 
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(3) The heat produced by vibration. 

(4) The elastic hysteresis loop. 

A complete physical theory would have to account for' 
all the observed results in all the preceding effects, For 
the purpose of flutter calculations, however, a complete 
physical theory is unnecessary". ' For "example , the exact 
shape of the hysteresis loop may he unaccounted for. It is 
sufficient to express the damping in terms of certain ef- 
fective parameters that can he conveniently incorporated 
into the flutter analysis. The present report reviews 
some of the results of— other investigators on the physical 
laws of damping but makes no attempt to find the true law 
of damping in the flutter models. ' 

The results of damping experiments are often expressed 
in terms of the loss of energy per cycle during "vibration. 
This energy can be considered as consisting of two parts: 
the part that is independent of frequency and the part 
that depends on frequency. The • independent part has been 
termed "statical hysteresis"; the dependent part, "hered- 
itary hysteresis." The results of numerous investigations 
indicate that, f-or most solids, the -statical hysteresis 
accounts for practically all of the internal damping. For 
example, Kimball and Lovell (reference 1) found that, in 
various metals, rubber, glass, celluloid, and maple wood, 
the loss of energy per cycle was independent of" frequency 
and could be represented with sufficient accuracy for most 
purposes by a constant times the square of the stress am- 
plitude. Keulegan (reference 2) found that the loss of 
energy per cycle in Armco iron was' the same whether deter- 
mined from a stntic-hysterresi s test- or from the rate of 
decay of vibrations. Other investigators (references 3, 4, 
and 5) have likewise found the loss per cycle to be inde- 
pendent of the frequency. 

Although the damping of materials is conveniently ex- ■ 
pressed in terms of the loss- of energy per cycle, it is 
difficult to write an analytical equation of motion that 
will represent the observed damping properties. Nearly 
all of the publi shed-analyti cal treatments express the 
damping as. a force proportional to the velocity. The con- 
stant of proportionality is thus a measure" of the equiva- 
lent viscous damping. This constant may then be regrad- 
ed as a function of amplitude and frequency. The use of 
the analysis of viscous damping for cases of other types 
of damping is based upon the approximation of assuming 
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that the hysteresis loop, can he replaced by an ellipse 
having the same ampli tude . and the. same area as the .actual 
loop. The damping coefficient determining this .equiva- 
lent ellipse is then used -to express the results of exper- 
iments. Ton Schlippe (reference 6) has applied this meth- 
od to the analysis of internal damping. 

Another analytical method of describing internal damp- 
ing has been used in reference 7. If the displacement is 
represented by a complex variable, a damping force propor- 
tional to amplitude but independent of frequency can be. 
represented by a complex stiffness constant. If this com- 
plex stiffness is written in the form k(l + ig), then g 
is a nondimensional damping coefficient. In the following 
analysis, the relation of the coefficient g to the vis- 
cous damping coefficient is shown and the principal formu- 
las for use in evaluating these coefficients from actual 
data are derived. 

The basis of the analysis of vibrations in continuous 
structures is the theory of normal modes. A normal mode 
means a type of vibration in which each particle of a struc- 
ture vibrates in simple harmonic motion with the same 
frequency and passes through its equilibrium position at 
the same time. The important property of normal modes that 
makes them useful in vibration is .the fact' that, when they 
exist, any possible type of vibration of a system can be 
represented by the superposition of vibrations in each of 
the normal modes and each normal mode can be treated inde- 
pendently as a system of one degree of freedom. It has 
been shown by various writers (see, for example, reference 
8, pp. 107-108) that normal modes will, certainly exist if 
there is no damping and if the potential and the kinetic 
energies are quadratic functions of the coordinates and 
the velocities of the system. In the theory of vibrations 
the amplitudes are usually assumed to he small enough that, 
in the expressions for the energies, all but the quadratic 
terms may be neglected. Even with damping, normal modes 
will still exist under certain conditions. It can also be 
shown that, when the damping is small, the theory of normal 
modes always gives a good approximation to the actual vi- 
brat i ons . 

When an external periodic force is applied to a struc- 
ture, all the normal modes are excited to a greater or a 
lesser extent. But when the amount of damping is small and 
the frequency of the applied force is near to one of the 
resonant frequencies of the structure, one of the normal 
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modes "becomes predominant in comparison with all of th«- 
others. Under this condition, the analysis for one decree 
of freedom will provide a good approximation for the re- 
sponse curve in that particular range of frequencies. The 



aon curve tor ine various modes are Known, inis innerier- 
nco effect can be evaluated. A number of methods are de- 
scribed for experimentally evaluating this effect. Simple 
theoretical expressions for the exact response curves have 
also been- found for the cases of a uniform cantilever beam 
in bending and in torsion. . : 

It is to be TLoted that measurements of the typo . con- 
sidered in this report yield thd total damping, includi-ng 
internal damping, air damping , losses in the supports, 
and whatever typos of damping are present. 

SYMBOLS 

m , mass. 
f(x)j restoring force. 

k, spring constant of a vibrating system. 
g(x), deviation of restoring force from Hooko's .law. 
F, F 0 , applied fxrrce';. instanianeous and maximum value. 

t , -time . 

x, displacement in system of one degree "of freedom, 

or position coordinate^ in a continuous structure. 

ATC , loss of energy or work dono -per cycle. 
A n • B n » °n » D n » constants. 

c, damping force per unit velocity. 
■ uu, angular frequency. 



u> 0 



, - natural angular frequency with no damping, ,J k /to . 
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A 
T 

5 > 5 n 
c cr 
a x 

M 
€ 
K 
H 

' u(x) 

In 
X n (x) 

0? 

V 

" * 

. a n > *>n 
Qn 

D m 
P 



natural angular frequency with -damping 



Ui, 



- X 2 . 



angular frequency of maximum response. 

c/2m. . • • 

period of free vibrations, 2tt/u 1 . 

nondi mensi onal damping parameter, c/mo> 0 . 

value of c for critical damping, 2miu 0 . 
1 x n 

— log ~ ; rrai is the logarithmic decrement. 
Trn in 

constant . 

damping coefficient used in reference 7. 
spring constant of a vibrator, 
radius of crank arm. 

amplitude corrected for effect of coupling with 
vibrator. 

displacement function for a continuous structure, 
generalized coordinate, 
normal function, 
kinetic .energy, 
potential energy, 
■dissipation function. 
c n , coefficients of mass, damping, and stiffness, 
generalized force. 

interval of distance along a beam, 
denominator of ■ mth term in equation (36). 
mass density. 
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polar moment of inertia of. section. 

torsion modulus of section, 

shear modulus of material. 

angular displacement. 

Young's modulus. 

moment of inertia of secti.on. 

area of sect-ion. 



(2n - l) ^/o)/o) n for large values of n. 



ANALYSIS 
Equivalent Viscous Damping 



The analysis of viscous damping is often applied to 
systems having a different physical law of damping. For- 
tunately, this convenient hut inexact analysis can give 
useful information "because certain approximations are jus- 
tified when the amount of damping, is small. The "basis of . 
this method follows. In viscous damping, the hysteresis 
loop is an ellipse. Corresponding to any - 'other type of 
damping, an equivalent viscous damping coefficient can he 
defined such that the hyst~e~resis ellipse will have the 
same amplitude and the same area as the actual Loop. The 
parameter characterizing this equivalent ^ellipse can then 
he used as a measure of damping, 

Consider a typical elastic hystere-sis loop (fig. l). 
This curve is seldom directly measured, "but a typical shape 
can he inferred for the purpose of- this discussion. Now 
the area of this loop is a measure of the" energy dissipat- 
ed per cycle, that is, a measure of. the damping. If the 
amount of damping is small, the loop must he narrow. Prom 
Hooke 1 s law, the mean slope should h-e approximately con- 
stant. The equation of motion can he written 

is + f(x) = F sin u)t (l) 
mx + g(x) + kx = T Q sin (jut 
f (x) = kx + g(x) 
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where ki is the elastic force corresponding to Hooke's 
law and g(x) is an undetermined function to take account 
of the damping. 



The energy dissipated per cycle is 

stt/uj 

dt 



AW = J" [kx + g(x)] If dt (2) 



and, whatever the actual law of damping may he, the ampli- 
tude of steady forced vibrations is determined by the con- 
dition that the work done per cycle by the external force 
is equal to the dissipation of energy by the damping. 

3tt/o) btt/o) 



* If d* = / f(x) |f dt = ATT (3) 

dt J dt 



Now the deflection will be very nearly a sine function of 
the time but, for generality, a Fourier series will be as- 
sumed. 

Let 

CO 

x = S A n sin nUJt + B n cos ncut 
n = i 

The phase can be adjusted arbitrarily to make B x - = 0. The 
velocity is - 

dx 00 

— -r = £ A n n(JU cos nuut - 3 n nuj sin ncut 
at n = i 

All of the coefficients except A x will be small when 
the amount of damping is small. Similarly, for a given am- 
plitude and a given frequency, the function g(x) can be 
expressed as a Fourier series. 



g(x) = 2 G sin nujt + D-n cos ncut 
n= l 11 " 



The elastic force . kx does not contribute to the damping 
and may be omitted from the expression for AW. 
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Then 




/27T/U) 
CO 
2 (C n sin nujt + IL cos nu>t ) 

CO 

X 2 (nujA n cos niwt - nU)B n sin nuut) dt 



= rr ( D a A x + 2D 2 A a - 2 C a - -K- . . . ) ( 4 ) 

All the cross products of sines. and cosines and the 
products involving two different values of n vanish in 
the integration. The only remaining terms are products of 
oach Fourier component of g(x) multiplied "by the corre- 
sponding; Fourier componentr—of the velocity. Thus, if the 
second harmonic and the higher terms in the Fourier series 
for x and g are small quantities of the first order, 
their contribution t-o the value of AW will involve only 
small quantities of the second order. For the definition 
of tho equivalent viscous damping, these higher order terms 
are neglect-ed in comparison with the term in Ai . 

Then 

AW = -it DiAx 

Under this assumption, the loss of energy per cycle depends 
only upon the Fourier component of— the restoring force in 
phase with the velocity. But it is just this component 
that corresponds to the term far viscous damping. The 
equivalent restoring force is 
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f (x) = kA x sin (jut + D^^ cos ait 

= kA x sin (jut + ~— cos out (6) 

This function gives a hysteresis loop in the form of an 
ellipse. The coefficient of viscous damping, c, which 
gives the same energy loss per cycle is given "by 

dx 

0 TZ = couAi COS U)t 



dt 



Hence 



AW 

— r— COS U)t 

-rrAi 



AW 



ttojA, 3 



(?) 



Derivations of formulas for One D e gree of Freedom 

In the following analysis, it lias been recognized that 

the quantity ui 0 = ^/Ic/m , the natural frequency without 
damping, is not an observable hut a conceptual quantity. 
Actual measurements can yield only the natural frequency 
with damping and the frequency of maximum response. Hence, 
the equations for determining' damping parameters should 
contain only the observable frequencies and not the fre- 
quency u) 0 . 

Fr ee vibration .- The equation of motion of a vibrating 
system with viscous damping will be taken in the form 

m'x + ci + ki = 0 ( 8 ) 



The solution is 



wher e 



-At . 

x = x Q e sin (juj.t 



/~~a ~a 
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The damping will "be expressed in terms of a nondimen si onal 
parameter defined "by 

This parameter is. simply related, to the fraction of crit- 
ical damping. The condition for critical damping is tUx =* 2: 
hence . . . . . 

„ 0 . X . |«. 

c _ c _6 
c cr 2bu) 0 "2 

The amplitude after n complete cycles of free ribration 
is ' 



x n = x o e - ~ A o 



n ctt 

-nXT " raSil 



1 , X Q 



= x Q e (10) 



mail Trn >c x n 



and since 



2 / „ 2 



5 1.x 



/ 



— 10?: (12) 



g- Trn ~ & x n 



Define 



Then 



trn e Xn 



a 2 
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Prom this equation it is apparent that, when 5 is small, 

8 ~ a x 

The quantity usually called the logarithmic decrement is 
equal to -na.^ in the present notation. 

For ged vibration produced b y a •periodic force of con- 
stant amplitu de. ~ In terms of the notation of complex var- 
iables, the equation of motion of forced vibration produced 
by a periodic force of constant amplitude is 

mz + cx + kx = P 0 e iUJt (13) 

The solution for a steady state is obtained by making the 
substitution 

icut 



Then 



Is. Za 

m m 



(14) 



Cw 0 s - O + i ^? (cu 0 a - a) 3 ) + icuu> 0 S 
or, for the absolute magnitude of x 0 , 



o 
m 



0 r— 1 — s 



(15) 



o / V i,,. 3 J ' u u; 0 2 



The frequency 0) m of the impressed force for which x D 
is a maximum is given by 

diu 2 IA w Q s ' tu 0 2 J uj 0 2 V cu 0 2 / tu 0 a 

i ^ a if 



(16) 
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The amplitude at resonance is 

J 



F o 



_o 

m . . 



2 



By usb. .of equations (15.) and (16), there is obtained 



a 

8 =2 



wher e 



(i - — 1 ) . . (18) 



2 >" 

V 1 U . B V 



a 2 

a 2 g- 



'5a; - i 



Forced vibratio n -pr oduced b,y a f orce p roportional to. 

yjl.*r The -equation of motion of forced vibration produced 
by a force proportional to uj 2 is 

mx + cx +- kx = MUJ 2 e itUt (19) 

where M is a constant. The amplitude of steady vi ora- 
tion is ?iven by 

M £2 Ma 

x 0 . ^- ■ 5.^ (ao) 



y - o'"'*^- j («>„" - »!■)"+»■ 



12 3 

m 



The frequency u> m of the impressed force for which x 0 
is a maximum is ^iyen by 
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duj* 



o - 



M 
m 



J~, 2 2 2 a 

J (UJ Q - w ) + 8 UJ % 



i - «?Y- 2(o> 0 2 - u>V + 5 2 ujo 2 
d. m r 



L' 



(UJ 0 S - U) S ) + 6 2 UU 0 a U) £ 



from which, is obtained 



UJ. 



m 



1 - 



3/S 



= 0 



(21) 



The amplitude at resonance is 



*m = 



m m 



•M' 

m 



y («.„» + 6 



22a 
^m 



■•A - ? 



(22) 



By use of equations (20) and (22), there is obtained 



(23) 



where 



a. 



2 _ 



Properties of Damping ' Co effi ci ent , g 

Theodorsen and fiarricfc (reference 7) hare used a com- 
plex stiffness function to describe the phenomena of inter- 
nal damping, The properties of this function are herein 
developed. 



The equation of motion of a vibrating system can be 
written in the form 
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nix + ik (1 + ig) = P 0 e iU)t (24) 



The use of the imaginary term ig implies that the dis- 
placement is a sine function of time and that there is a 
component of force proportional to the amplitude and inde- 
pendent of frequency but in phase with the velocity. The 
solution of equation (24) is 



x = 



P Hot 
k 6 



i - J»l + _ u 



(25) 



This equation agrees with Schlippe's analysis (reference 
6) and is similar to the result just given for viscous 
damping except that the damping term is ig instead of 
i((jj/uu 0 )5. Thus, the numerical value of g is nearly equal 
to the value of- 8 obtained by measurements close to a 
resonant peak. 

The frequency of maximum response is obtained as be- 
fore by differentiating x with respect to U) 2 . The re- 
sult for this case is 

OJ _ tu_ . 
m o 

Hence • 



I 
k 

x 0 = sr 



1 - + ig 



k_ 



(l - -%) 

g 3 „ A J*3Ll (26) 

(fsY - i 

\x 0 / 



This equation shows that g is given by-" the same expres- 
sion as the quanti-ty a 2 of equation (18). 
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Further Details of Practical Importance 

In the elementary analysis, the applied force is usu- 
ally assumed to be independent of the amplitudes produced 
in the vibrating system. ; in actual apparatus, however, 
the motion of the point where the external force is applied 
frequently affects the value of the force, so that the ac- 
tual force transmitted to the structure changes with the 
amplitude of vibration. When, the force is produced by a 
spring of stiffness E fastened to a crank of radius R 
(see fig. 2(a)), the force .on the structure is 
- K(x - R sin tut) and the equation of motion is: 



hi 



+ kx = - K(x - R sin wt ) (27) 



Then 



1 -a) 2 m + k 1 



x KR sin tut R sin tut 



1 - 1 + i 



x Q x Q ' R 



- 1 _ _1 1 

x 0 1 x 0 R 



(28) 



Similarly, for a rotating mass fastened to the structure 
(see fig. 2(b)), the equation of motion is: 

mx + kx = - M -~~e (x + R sin out) (29) 
at 

The solution is : 

1 -mou 2 + k 1 



Mu) s R sin cut R sin out 



-o -*-o 



1 
R 



1 1.1 



c 0 ' x Q R 



+ ± (30) 



The corrected amplitude, x Q ' , corresponds to the ampli- 
tudes given by theoretical equations in other sections of 
this report. Equations (28) and (30) thus provide a method 
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of correcting the isa-sured amplitudes for the effect of a 
vi orator having a small stroke. The natural frequencies 
are changed to -- . ; 



/ k + K 
v m 

and , 

jmz ■ 

V m + M 

for the two types of vibrators shown in figure 2(a) and 
2(h), respectively. 



Derivation of Formulas for Continuous Structures 



General case." The general theory of. vibration has 
been presented "by a number of authors » (See, fox example, 
reference 8, chs. IT and V.) This theory deals with the 
problem of finding the normal coordinates of a system, 
that is, the coordinates in terms of which the equations 
of motion have only one coordinate occurring in each equa- 
tion. The existence of such coordinates for systems with 
damping has been discussed by Lord Rayleigh, who shows 
that these coordinates exist for certain distributions of~ 
damping and that, in all cases of small damping, the errors 
introduced by assuming normal modes are of the second order. 

A brief outline of the m-ethod of normal coordinates 
is given. The application of this method to the determi- 
nation of damping in each mode of a cantilever beam is then 
discussed. 

Suppose that the functions X^x), giving the shape 
of the deflection curve for each-mode , are known. Then 
the displacements corresponding to any motion of the syst-em 
can be expressed as a summation of the displacements in 
each normal mode 



u(x) = E q n X n (x) (31) 
n 

According to Lord Rayleigh (reference 8, pp. 130-131), 
in damped systems for which normal modes exist, the kinet- 
ic and the potential energy functions, T .and V, and the 
dissipation function, F, can be expressed as sums of 
squares of generalized coordinates or velocities. 



Then 
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21 =» &1 4i S + aa q^ + ... + a n q n a + . . . 

2F = fcj. q.i 2 + "ba^a 2 + ••• + ^n^n 2 + ••• 
2 V = c 1 q 1 a + c 2 q 2 2 + ... + c n q^ s + ... 

The equations of motion will "be obtained "by the use of 
La?ran?e'e equations. 

dt ^3q n ' oq r _ 3q n 

a n*n + + °n*n = ^n < 33) 

where Q, n is the generalized force corresponding to q n . 
The form of Q n is found from the relation 



Qn = 



8W 



n " Sq n 

where 8W is the work done on the system hy the external 
force during a displacement 6q n - 

For example, in the- case of a force IP = Pge^ ap- 
plied at the point x = x x , . 

5¥ = %l 5 <ln = J8v - = ^o 6 ^ 5( ln 

Eence 

Q n = X^xJ F 0 e iujt (34) 
The solution of the equation .of motion is 

_ X n (x 1 ) F Q e^ b 
q^ — g 

-Lt) a n + iuTo n + c n 



X n (x 1 ) F Q e iUJt 



a^U) 2 



^on ^on 7 



(35) 
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where 

Won - ari 



!he displacements are .then given hy 



«(«) - E 4n X n . S *» ( *' ^ ( V> (36) 



In the summation, it is necessary to sum the real and the 
imaginary parts separately. If only real quantities are 
used, the equations "become 

Q n = X n (xJ F 0 sin out 

q n = A n sin u>t + B n cos tut 

X n (x x ) Fo (l - /a) „ 



A 




n 2 3 

an ^ K 1 - + ^ 6n ) J 

■a(x) = £ A n X n sin uit + B^. X n cos ujt (37) 
n 

Equation (36) is the connecting link Jiotiveon the analy- 
sis of vibration in one decree of freedom and in continuous 
structures. Each term of- this equation has the same form 
as the solution for a system having one degree of freedom. 
When a force is applied to a structure, all the mod e s are 
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excited to -a greater or a lesser extent and the problem is 
to find the amounts of damping associated with each mode. 
There are several ways of solving this problem. 

For any case in which the normal functions are known, 
the damping in different modes can be separated by an in- 
tegration. If both sides of equation (36) are multiplied 
by X m (x) and integrated over the length, then all but 
one term in the summation will vanish. 

u(x) = 2 q n X n 
n 

/I pi 
X m (x) u(x) dx = / X m (x) S q n X n (x) dx 

o o 

= <L m / x m S (x) dx 



/ l ^,(x) u(x) dx ^(x,) F 0 e iujt 

^m = — n : = i — > - a ~ : < 38 > 

f l X m 2 (x) dx a m % m (l - ^- + i ^ 8 m ) 

o \ ^om ^om ✓ 



The- quantity q m obtained by graphical or numerical inte- 
gration of the measured amplitudes u(x) can then be used 
as though it were the amplitude in a system having only 
one degree of freedom. The process of evaluating this in- 
tegral in practical cases consists in measuring the ampli- 
tude at a number of points and evaluating the sum: 

S X^Cxi ) u(x ± ) AX£ 

q m = (39) 

Z X m (x ± ) Axi 
x 

In the present tests, the procedure was simplified 
still further. The disturbance from two of the modes can 
be eliminated by applying the force .at a node of one of 
them and measuring amplitudes at a node of the othe.r. The 
justification of this method follows immediately from equa- 
tion (36), For measurements in the range of frequencies 
near to UJo m . this equation can be written: 
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u( t ) - ^-i** "* * n -.i<*i> F o eiUit + X m (x 3 ) X B ( Xl ) F 0 e luA 

+ ^_^_^±^2_^ + eDall term8 (40) 

where D n is written for the denominator in the rath terra. 
If Zi and x 2 are chosen such that 

then the .jaiapli tude is given by 

u(x 2 )=0+ — : — a ~ : + Of snail t"erns (41) 



an w on (1 - + i — S n ) 

> or.i on * 



on on 

Another method of finding the damping is by evaluat- 
ing the infinite series given "by equation (36) to find the 
resultant response curve*. A comparison of this exact 
curve with the response ... curve for one dsCr.e_e' of freedom 
shows how to correct the analysis for one degree of free- 
dom to take account of all the disturbing modes. Expres- 
sions for the sum of the infinite series (equation (36)) 
have been found for the cases of torsion and "bending of a 
cantilever beam with the vibrator at the tip. 

To rsion- of a. uniform bea m.- Consider a cantilever beam 
excited in torsional vibrations by an oscillating torque 
applied at the fre-e end. 
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■Hence 

,1 



a n = J ■ £ P I p. x n 2 ^ x) dx 



o 

b n = 0 



\ ( M 



3 



Suppose that the applied torque is supplied by a rotating 
unbalanced mass at the tip of the beam. Then the torque 
i s 

P = Mu> 2 sin uut 
The generalized force is 

Q n = X n (l) Mcu 2 sin cut (43) 

Then 

^(l) Mw 2 sin u>t 

6 U) = S q n X n = 2 _ (44) 

n - 



on 

The normal functions are known for the case of constant Ip, 
X n U) = sin (2n - l) ||; n = 1, 2, 3, ... (45) 

9(1) = g. sin S 



UJ 

Put 

_£n. - 2n - 1 



8(1) = }Ll^J>* 2 1 _ ( 46 ) 

|P' lI p tt f (2n - l) 3 cu x a _^ x "j 
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The sum of this series can "be evaluated by the following 
method. Using the relation (reference 9) 



tan x = - 2 — -~ 

-co x - (n + § ) TT 



(47) 



there is obtained 

\ 



x tan x = 

l 



(n - 



i) J - 1 (n. - * 5 + 1 



CO 

i 



(n 



, v" IT 

f) - 
X 



(48) 



- 1 



Put ix = ju 

TT Ulj. 

Then 

CO 



T JUK TTOJ _ v -- 

^ tan ^— _ E (2n _ 1)2 a 

1 

oi 2 



q / 7 \ M TTO) , 10TT 

9 0 (-W = -i -- tan — - (49) 

0 i p I I p 4u)i ujj.2 



In figure 3 the value of 

e 0 (i) s p i i. 



has been plotted against uj/iUj. fxrr the first two peaks. 
In figure- 4-5 the exact response curve is compared with the 
curve for one degree of freedom. 

Bending of a uniform beam.- By a method similar to 
that used for torsion, the following aquations are obtained, 
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T = 



p A u 2 dx = E. q. 



n 



g p A X^ dx 



dx 



(50) 



a = J* 4 p AX n 2 dx 
•b = 0 



c = 



Suppose that the applied force is produced "by a ro- 
tating unbalanced mass at the tip of the "beam. Then 



Q n = Xj^l) Muj 2 sin tut 



u(x) = S ci 2L = E — 



X n Gc) X a (l ) Mou 2 sin cut 



u(x) = E 



n n n n a r , ~ . x / ou^n 

M X n (x) sin tut 



Put 



y = E 



and define a new function 
oo ' , 



y' - E — 



n=i (2n - l) 



(51) 



(52) 



(53) 



(54) 
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yi = (tan Sg! - tanh Sgl) (55) 

If cp is defined "by (2n - 1 ) "/qp" = w n /«J for large 
values of n, the function y' is practically identical 
with the function y, except f-or the fir-gt two terms of 
the series. By use of this function y' and suitable cor- 
rections fD-r the first two terms of the serie-s, the re- 
sponse curve has been calculated for the fij^st three "bend- 
in? modes of a beam and is shown in figure 5. Figure 6 
shows a comparison of these exact resonance peaks with the 
corresponding curve for one degree, of freedom. 



APPARATUS 



The apparatus used for these teats was very simple. 
The applied force was produced in part of the tests by a 
small rotating unbalanced mass, made by tap ping a screw in- 
to the side of- a rotating shaft, fitted in iight bearings 
and clamped to the structure that- was being te-stod. The 
shaft was turned by a small electric motor fitted with & 
fl exible coupling and suitable ^ears; the frequencies were 
measured by an electric tachometer. For the rest of the 
tests, the rotating mass was replaced by a small crank 
coupled to the structure by a rubber band. 

The amplitudes were measured by observing the magni- 
fied image of the filament of a small bulb produced by a 
lens of 1/4-inch diameter and 2-^3/4-inch focal length- held 
on the structure by a small brass mounting. 
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The models tested were, the wings used in. the flutter 
investigation, (reference 7) and were in the form of can- 
tilever teams 6 feet 9 inches long with symmetrical- 
airfoil cross sections. Model 1 was rectangular with 12- 
inch chord and l/2-inch maximum thickness. It was made of 
duralumin with closely drilled l/2-inch holes- and was cov- 
ered by a 0,006-inch sheet of- duralumin. Models 2A , 20, 
3, and 4 were made of solid duralumin with two rows of 
chordwise slits to decrease the torsional stiffness. Mod- 
els 6 and 7 were made "by covering a balsa structure with 
l/l6-inch mahogany. 

EXPERIMENTAL- RESULTS 



Three methods of measuring damping were used in these 
tests. The damping in the first bending mode of each 
model was determined by measuring .the rate of decay of 
free vibration. For the higher modes, some of the tests 
were made with the rotating-mass type and some with the 
crank-and-rubber-band type of vibrator. In order to min- 
imize the influence of normal modes other than the one 
being investigated, the force was applied and the ampli- 
tude was measured at the most suitable positions. lor ex- 
ample, in the measurement of the damping in the third 
bending mode, the nearest ' di sturbing frequencies are the 
torsion and the fourth bending. In this case, the force 
was applied to the tip of the beam at the position of the 
node in torsion and the amplitudes were measured at the 
edge of the beam at the position of- a node of the fourth 
bending mode. 

Figure 7 shows plots of the response curves of the 
models tested. In these tests, it was thought unnecessary 
to try to obtain the complete functional dependence of 
damping on amplitude and other variables. Consequently, 
the response curves were analyzed only to the extent of 
finding a representative value of the damping parameter 8. 
It will be noted that the results for all models of a given 
type of construction gave values of 6 within a range of 
about a factor of 2. In table I are given the numerical 
results for 8. These values were computed from the data 
of figure 7 by use of equations (12), (18), and (23). The 
value of the amplitude at a frequency ratio of 1.1 was 
taken for determining the best single representative value 
of 8. 
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conclusions 



1. In the determination of the damping of structures 
"by means- of the" shape of the response curve obtained hy 
applying an alternating load at one point of the .structure , 
the use of the analysis- for one degree of freedom is jus- 
tified when the following conditions are met: 

(a) The damping- is small. 

(b) The points of applying the force and measuring 

the amplitudes are appropriate from consider- 
ations of disturhing normal modes. 

(c) Only amplitudes close to a resonant peak are 

used to- determine the nondimen si onal dampi-ng 
parameter, 8. 

2. When the normal functions for a structure are 
known, the damping in the different modes can "be separate- 
ly determined from the measured amplitudes at several 

.points along the structure. 

?■ . The -measured values of 6 for a homogeneous 
structure such as a cantilever beam of -duralumin are ap- 
proximately equal in the different modes of vibration. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, 7a., December 21, 1939. 
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TA3LS I 1 

Tallies of the Damping Parameter, 8 







Mode 






Model 


3endin? 


Torsion 






1 


2 . . 




. 4 


1 


1 .. 


jj ll r a ±-uin l n 
skin 


0.062 


0.083 


0.1130 




0.087 


2A 


Solid 
duralumin 


. 0050 


.014 






.0087 


20 


do . 


. 0050 


. .0056 


.0073 ■ 




."OIB 


3 


do . 


. 0045 


.0075 


.0040 ■ 


6.0047 


— — * 


4 


do . 


. 0089 


. 0080 


. 0050 


' .0053 


.0034 


6 . 


Wood 


. 020 




. 0174 


. 030 


.028 


7 


do . 


. 021 








. 020 
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Figure 1,- A typical elastic hysteresis loop. 




Figure 2.- Simple mechanical circuits to illustrate 
the effect of a vibrator. 
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Figure 6a to o.- Comparison of response curve in bending with curve for one degree of freedom. 
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Figure 6.- Concluded. 
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Fig. 7 



6=0.014 



6=0.0087 



8=0.0070 



(a) 



(b) 



Co). 



6=0.0070 



6=0.015 



6= 



0.008 



(f) 



GO 



(e) 




8=0.0050 




6=0.0055 



8=0.0034 



(8) 



(h) 



(i) 



.8 



1.0 1.1 1.2 .8 



1.0 1.1 1.2 



.9 1.0 1.1 1.2 .8 

(a)Model 2A; second bending (b)Model 2A;torsion mode; (c)Model 3; second bending mode 

rotating-weight method. 



mode; rotating-weight 
method. 



rotating-weight method. 



(d)Model 3; third bending (e)Model 3;fourth bending (f )Model 4; Beoond bending mode 



mode; rotating-weight 
method 



mode; rotating-weight 
method. 



rubber-band method. 



(g)Model 4; third bending (h)Model 4; fourth bending (i)Model 4; torsion mode; 
mode; rubber-band method. mode; rubber-band method. rubber-band method. 



Figure 7, a-i. -Typical experimental response curves. 



